function irf = RPfigures_b(pars,T)
%Compute IRF woth an ARMA(2,1) structure
%Read ARMA parameters
eta = pars(1); rho = pars(2); mu = pars(3); nu = pars(4);
%Compute IRFs
irf(1) = mu; irf(2) = (eta + rho)*irf(1) + nu;
for i=3:T; irf(i) = (eta + rho)*irf(i - 1) - eta*rho*irf(i-2); end
